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The partner mode with respect to a vacuum state for a given mode (like that corresponding to one 
of the thermal particles emitted by a black hole) is defined and calculated. The partner modes are 
explicitly calculated for a number of cases, in particular for the modes corresponding to a particle 
detector being excited by turn-on/turn-off transients, or with the thermal particles emitted by the 
accelerated mirror model for black hole evaporation. One of the key results is that the partner mode 
in general is just a vacuum fluctuation, and one can have the partner mode be located in a region 
where the state cannot be distinguished from the vacuum state by any series of local measurements, 
including the energy density. I.e., ’’information” (the correlations with the thermal emissions) need 
not be associated with any energy transport. The idea that black holes emit huge amounts of energy 
in their last stages because of all the information which must be emitted under the assumption of 
black-hole unitarity is found not necessarily to be the case. 

PACS numbers: 04.70.Dy, 04.62.+v, 04.60.-m. 


I. INTRODUCTION 

Quantum mechanics around a black hole has been one 
of the most exciting and puzzling aspects of theoreti¬ 
cal physics in the past half century [lj. One of the is¬ 
sues has been that of’’information” and how information 
is carried. Black-hole unitarity [1] is the belief that in 
the space-time outside the black hole, the evolution of a 
quantum field from a time before the black hole formed to 
after the black hole completely evaporated must be uni¬ 
tary - initial states of the quantum field map uniquely 
to final states of the quantum field. If one believes in 
” black hole unitarity” (and in this paper we are agnostic 
about that belief) then there must exist correlations be¬ 
tween the early Hawking evaporation emission from the 
black hole, and late time emission. For any particle emit¬ 
ted early on, some correlation between this early mission 
and the field later on must exist. Given a mode which 
carries away thermal particles in the early stages, there 
must be ’’partner modes” which occur later which are 
correlated with these early modes in order that ’’unitar¬ 
ity” be preserved. 

The characterization of these partners thus becomes 
important. In section eh we define the partner mode, 
uniquely in some cases of pure two-mode squeezing. In 
section mu we show that any particle detection measure¬ 
ment of the field also has a partner, even in the case 
where the detector is stationary but is switched on and 
off. In section [V] we look at the partner in the case of the 
accelerated mirror model of black hole thermal emission. 
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One of the surprising results is that the partner need 
not be located near the original mode, but can be located 
in distant regions of the space-time. While recognized in 
the correlations between the field inside and outside the 
black hole in the Hawking evaporation process 3, this 
is a general feature of the partner modes. Our results 
are somewhat related to to recent observations that the 
long range entanglement in the vacuum can be used to 
entangle other systems even in spatially separated regions 
0 and to energy teleportation studies @ ■ 

First, let us specify how the partner mode can be de¬ 
fined. To this end, we demand the two conditions: 

A The reduced density matrix of the Hawking plus 
partner modes obtained by integrating out all other de¬ 
grees of freedom should be a pure state. Since the total 
state (the initial vacuum) is a pure state, this is equiva¬ 
lent to vanishing entanglement between the Hawking plus 
partner mode on the one hand and the rest of the system 
on the other hand. 

However, this requirement alone does not define the 
partner mode uniquely (see below). For example, one 
could envisage a single-mode squeezing operation and 
phase rotation acting on the partner mode, which does 
not change the purity of the combined state (Hawking 
plus partner). Specifying the partner mode uniquely re¬ 
quires a second condition. There are several reasonable 
options, here we list some possibilities: 

B1 The quantum state after absorbing (annihilating) 
one partner particle should be (up to normalization due 
to possibly different probabilities) the same state as after 
creating one Hawking particle. This corresponds to the 
intuitive picture that the Hawking and partner particles 
always come in pairs. 

Alternatively, one could implement the idea that 
Hawking and partner particles always come in pairs by 
imposing the requirement the other way around: 
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B2 The quantum state after absorbing one Hawking 
particle should be (again up to normalization) the same 
state as after creating one partner particle. 

As we shall see below, condition B1 can always be sat¬ 
isfied - unless the Hawking mode contains single-mode 
squeezing only and thus there would be no need for a 
partner particle at all - whereas the requirement B2 can 
only be fulfilled if the single-mode squeezing of the Hawk¬ 
ing mode is small enough. As another option, we could 
demand that the probabilities for detecting Hawking and 
partner particles should be the same - treating these two 
modes on a symmetric footing. As it turns out, in the sce¬ 
narios we are interested in below (pure two-mode squeez¬ 
ing) , all these requirements yield the same answer for the 
partner particle. 


II. DEFINITION OF PARTNER PARTICLE 


Now let us show how to satisfy these requirements. As 
a most general ansatz, we decompose the Hawking mode 

a# = J dk (a* k a k + , (1) 

into creation and annihilation operators a' k and d k 


[dk,dk>] = 




= 0, 


dk,d\, 


= S(k,k'), (2) 


where |ny| 2 = {ny|ray} = 1 = |nj_| 2 and ny _L n±, i.e., 
{«,|||nj_} = 0. If /3_l = 0, we would have pure single¬ 
mode squeezing and the Hawking mode itself would be 
in a pure state, i.e., there would be no need for a partner 
particle. In the general case f3± ^ 0, we can restrict 
ourselves to the two modes 

a|| = {n|||d} , a_L = {n_L|a} , (9) 

which satisfy the usual commutation relations due to 
Eq. |6]). These operators annihilate the initial vacuum 

d\\ |0) = fij_ |0) = 0, (10) 

and thus the reduced density matrix of these two modes 
is a pure state. Now the idea is that everything involving 
the Hawking mode dp and its partner mode ap will occur 
in the two-mode space spanned by dy and dp and their 
adjoints ajj and dp. As we show the the Appendix, this is 
actually the only way to satisfy requirement A. In terms 
of these operators, the Hawking mode is given by 

o-h = a*d|| + /3||dJ| + P±_a\ . (11) 

From [dp, dp] = 1 follows \a\ 2 — |/3||| 2 — |/3_l| 2 = 1. Note 
that one can make the Bogoliubov coefficients a, /3», and 
(3± real by absorbing their phases into dy, dj_, and dp. 

Following our strategy, we can make the following gen¬ 
eral ansatz for the partner particle 


defined with respect to the initial vacuum state 


ap = 7ydy+7*dj_ +(5||aJ| +dj_a 1 L . (12) 


Vfc dfc |0) = 0, (3) 

where k denotes some quantum number. 

For convenience, let us introduce the usual complex 
scalar product of two functions or vectors x and C via 

{X\C} = JdkxUk- (4) 

Accordingly, we define the projection of the initial anni¬ 
hilation operators a k onto one mode x via 


a x = {xl«} = J dk x* k a k , 

which gives the commutation relations 


*f] = 




= 0, 


a x ,d^ 


= {xlC} 


In this notation, the Hawking mode is given by 
dp = {a|a} + ({/3|a}) t = { a\a} + {a\/3} . 


(5) 


( 6 ) 


(7) 


Now let us introduce an orthonormal basis ny and n± in 
the subspace spanned by the two vectors a. and (3 


( 8 ) 


In this way, requirement A is automatically satisfied. 
Since dp should obey the usual commutation relation 
[dp, dp] = 1, the above Bogoliubov coefficients should 
satisfy | 7 y| 2 + |7_l| 2 — |d||| 2 — |Al| 2 = 1- Furthermore, 
since we want the two modes dp and dp to be in¬ 
dependent, i.e., [ dp,a} H \ = 0 = [dp, dp] as well as 
[dp, dp] = 0, we get the conditions 7 |*a = /3y<5y + (3]_S± 
and 7y/3y + 7 jL/?-L = a*<Sy. 

As mentioned above, these three equations do not spec¬ 
ify the four Bogoliubov coefficients for dp uniquely. We 
could still apply a single-mode squeezing/phase trans¬ 
formation dp —>• e lv cosh£ dp + e 1 ^ sinh£ dp within the 
dp-mode for arbitrary (real) values of <p, £, and ip with¬ 
out violating any of the conditions above. In order to fix 
this remaining degree of freedom, a second requirement 
is necessary - here, we discuss B1 and B2. 

Option B1 corresponds to choosing dp |0) oc a) H |0), 
i.e., 5||o:, which means <5j_ = 0 and dy = S. This then 
gives 7 y = (dyd/a and 7 p = (/3J 1 + /3]_)5/a such that the 
remaining Bogoliubov coefficient S can be determined by 
the unitarity condition |'y| 2 — |d| 2 = 1 up to a global 
phase. Writing this Bl-condition dp |0) oc dp |0) in the 
form dp |0) = r] dp |0) with some constant rj, we have 

(dp - rydp^ |0) = 0. 


a = any , /3 = j3\\n\\ + fi±n± , 


(13) 
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Thus, this linear combination ap — rja^ is composed of 
initial annihilation operators only. 

The other option B2 corresponds to ap |0) oc dp |0), 
i.e., -y||/3, which allows us to determine the Bogoliubov 
coefficients in a completely analogous manner. Note, 
however, that there is an important difference: As 
shown above, condition B1 can always be fulfilled unless 
f3± = 0, in which case we would have pure single-mode 
squeezing within the a//-mode and there would be no 
need for a partner mode. In contrast, requirement B2 
cannot be satisfied if the amount of single-mode squeez¬ 
ing becomes too large. 

In case of vanishing single-mode squeezing /3|| = 0, 
both requirements give the same partner mode 

ap = a*a± +/3al , (14) 


where a and /? can be made real because their phases 
can be absorbed into the definition of a\\ and a±. In this 
case, the initial vacuum state restricted to the two modes 
an and ap is a pure two-mode squeezed state 

|0) = exp fip&p — h.c. j |0)p F , (15) 

with respect to the zero-particle state \0) HP which is an¬ 
nihilated by ap and ap 


a h |0 ) HP — ap |0)pp — 0 , (16) 


where the squeezing parameter £ satisfies a = cosh£ 
and j3 = sinh£. As a result, the initial vacuum |0) can 
be viewed as a state containing pairs of particles in the 
modes ap and dp. After tracing out (averaging over) the 
partner mode, this squeezed state (1151) yields a thermal 
type density matrix for the Hawking mode 


1 

PH = ~z exp 


CLpCLp ( 

T H /’ 


(17) 


with the normalization Z ensuring Tr{/3p} = 1 and 

T “ = 2ii,(cLiie) • (18) 

which can be regarded as a dimensionless Hawking tem¬ 
perature. Due to the symmetric nature of the squeezed 
state (usd, the same applies to the reduced state of the 
partner particles (after tracing out the Hawking mode). 

Note that the mapping from the Hawking mode ap to 
its partner mode dp is not linear in general if ap has 
the partner mode dp and a' H has the partner mode dp, 
then the partner mode for /idp + va' Hl for example, is 
almost never /aap + va ' P , even if we have no single-mode 
squeezing in both cases. 


III. PARTNERS AND DETECTORS 

The idea of a partner particle has a broader applicabil¬ 
ity than just Hawking or acceleration (Unruh) radiation. 


Consider a model particle detector as suggested by Unruh 
Q and developed by De Witt Q. The detector is taken 
as occupying a single point in space-time with an inter¬ 
nal degree of freedom, often taken to be a spin degree, 
but could equally and more simply be taken to be a har¬ 
monic oscillator degree of freedom. The energy difference 
(in the rest frame of the detector) between the ground 
state and the first excited state is E. This is coupled to 
the quantum field of interest. This detector responds to 
specific degrees of freedom of the field, changing its state 
from ground to excited state, which is regarded as a de¬ 
tection. (I.e., if the detector is discovered at some time 
to be in its excited state, it must have absorbed energy 
and a particle from the field.) 

The interaction Lagrangian is given by (H = c = 1) 

Tint = e(r) q{r) d T §[t(r), x(t)\ , (19) 

where e(r) is the possibly time dependent coupling, 
f(r), x(t) is the trajectory of the detector in terms of the 
proper time along the path r. After quantization, the 
internal degree of freedom of the detector corresponds to 
the operator 


q(r) = ce~ iEr + c ! e iEr , (20) 

where c is the annihilation operator taking the detector 
from the first excited state of energy E to the ground 
state of zero energy. Note that if the det ecto r is a har¬ 
monic oscillator, then c could just be y/2E times the 
usual oscillator annihilation operator. The normaliza¬ 
tion of c is not important because it will cancel out in 
the following anyway. 

One can define a field operator associated with the 
detector by 


ap = Af 


dr e(r) e lEr d T $[t{r), x(r)], 


( 21 ) 


where M is chosen so as to make 


dp,a ] D 


= 1 . 


( 22 ) 


Furthermore, one can define a mode function associated 
with this operator by 






(23) 


Assuming that the field is (initially) in a vacuum state, 
the excitation probability of the detector will be given by 

( 0 | d^dp | 0 ) . . 

M 2 K ’ 

Hence ap corresponds to the Hawking mode ap. This 
mode, (j)p, is the mode that the detector absorbs when 
it is excited. This mode, by the above argument, has a 
partner mode, (j>p which is orthogonal to (f>p but is per¬ 
fectly entangled with cf>p in the vacuum state |0). If one 









4 


has a detector to measure various attributes of that part¬ 
ner mode, one will get vacuum values if one ignores the 
outcome of the measurements of the detector mode (f>p■ 
Since ap is a mixture of positive and negative (pseudo) 
norm vacuum modes, one will find a non-zero probability 
of finding it in the vacuum. But that probability would 
be the same as the probability if one looked into the 
vacuum without measurement of the Hawking/detector 
mode 4>d- There would of course be correlations between 
the ap and dp modes. I.e., if the Hawking mode were 
detected (the detector was found in its excited state), the 
detector measuring the partner would also be excited. 

Inserting the usual representation of a massless scalar 
field in 1+1 dimensional flat space-time 


<l(f,x) = 
we have 


dk 


d k e~' i (l fc l t_fex ) -(- e i{\k\t-kx) 


,(25) 



dr e(r) 


dk 


I\k\ / dt 

V 47 t \dr 


k dx\ 
\k\ dr ) 



-i(|fc|[t(r)-x(r)]+Br) 


- ai e 


t +(|fc|[t(T)-x(T)]-.E-r) 


)' 


(26) 


This allows us to read off the Bogoliubov coefficients of 
section eh expressing dp, the equivalent of the Hawking 
mode annihilation operator, in terms of the operators dk 
and d\ 


JV 


a h = 


dr e(r)\ 


e -i(|fc|[t(T)-x(r)]+Er) 


fa = 


0 i(|fc|[t(r)-x(r)]-£;r) 


k dx\ 


\k\ dr J 

(27) 

k dx\ 


\k\ dr J 

(28) 


In general, both the a coefficients and /3 coefficients will 
be non-zero, and a k will not be proportional to fa. Thus 
there is a partner mode. 

Let us now restrict attention to the case where the 
detector is at rest at x = 0, but e(t) is non-trivial. We 
choose e(t) such that 


Now, the term in the first line is just the positive fre¬ 
quency part of — fa D , which we can write as 


dk 


■ fae 


—i\k\t _ 


= - dt 


, </>p(t',x = 0) 

2iri(t — t' — i0+) 


J \Z4n\k\ 

Hence we can also write the partner mode as 

, (/>* D {t',x = 0) 


■ (32) 


<j)p(t,x = 0) = — cothr J dt' 


- tanh j 


dt 


2ni(t — t 1 — *0+) 
, fa D (t’,x = 0) 
2iri(t — t' + *0+) 


.(33) 


Thus in general (f>p(t, 0) will have a long tail falling off 
as 1/t for large |t|. Only if the moments J dtt n <j)p(t) 
are zero for all n will the partner fall off faster than any 
power. (If those are zero for all n < IV but non-zero 
thereafter, then </>p will fall off as l/|t| JV+1 .) 

Let us now give an example. Let us assume that 


e(t) = e 0 (i) + e Q (i - T)Xcos[3E(t - T)} , (34) 

where A is very small. Again E is the energy difference 
between the two states of the detector. The remaining 
function eo (t) is supposed to be a smooth switching func¬ 
tion. Furthermore let us assume that the Fourier trans¬ 
form of eo(t), namely ?o(w) is real and non-zero only in 
a compact region —E/4 < uj < E/4. The Fourier trans¬ 
form of e(t)e lEt , which occurs in the expression for ap, 
will have three peaks, one small one centered at —2 E, 
one large one at E and another small one at 4 E. Thus 
fa will have two small peaks with amplitude proportional 
to A at k ~ ±2 E, while a k will have two large peaks of 
amplitude 0(1) at k « ±E and two of amplitude A at 
k « +3.E. Because of the limited width of each of these 
peaks, none overlap, and (3 will be orthogonal to a. Thus 
the “detector mode” will be a pure two mode squeezed 
state. The /3-coefficient and thus tanhr will be of order 
A while cothr will be of order 1/A. The partner mode 
will have a temporal Fourier transform with a single peak 
centered at 2 E of amplitude 0(1), two smaller peaks of 
amplitude A at —E, and one of amplitude 0(A 2 ) at —3 E. 
Thus the partner mode will be approximately given by 


{a\(3}= J dka* k fa =0, (29) 

so that the detector mode corresponds to pure two-mode 
squeezing. Defining 


</> P (t,x = 0) = eo(t-T)e 2iEit - T) xO(l) 
+e 0 (t)e~ iEt x 0(A) 

+e 0 (t - T)e~ 3iE ^- T) x 0(A 2 ). (35) 


cosh 2 r = {a|o;} = J dk\a k \ 2 , 

sinh 2 r = {/3|/3} = J dk \fa\ 2 , (30) 

the partner mode is 

c/p = f —( f 3 ke - i ^ t - kx faothr+ 

J y/4'K\k\ A 

+a* k e i{Wt ~ kx) tanhr) . (31) 


The envelope of the partner mode will thus be dominated 
by eo (t — T), i.e., displaced from the detector mode by a 
time T. As a result, the partner mode will be centered 
around a time arbitrarily displaced from the maximum 
of the detector mode. Of course, this is somewhat mis¬ 
leading since the part oc eo (t — T) of detector mode which 
leads to detection (the fa part of the detector mode) and 
the a* part of the partner do overlap. 

However, if one chooses some other mode, (j>x which 
has an overlap with the partner but, let us assume, none 
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with the detector mode, there will be correlations be¬ 
tween measurements made on this mode and the out¬ 
comes of the detector measurements. 

We note that this example shows that partner modes 
are not a unique feature of black holes, or accelerated 
detectors. All detectors, which have a finite probability 
of detecting something in the state of interest, even if 
due to “switch on/off” transients, will have both a “de¬ 
tector mode” and a partner mode associated with them. 
If the partner is well separated from the Hawking mode 
(which it is if we are interested in the detection of radi¬ 
ation from say a black hole, where the partner is behind 
the horizon, and the other is far from the black hole), 
then any measurements made on the partner mode will 
give results indistinguishable from the results in that vac¬ 
uum state. There will however be correlations between 
the results for measurements on the partner and on the 
“Hawking mode” and not with any other modes orthog¬ 
onal to these two. If one were able to communicate be¬ 
tween the detectors, one could for example measure the 
partner whenever the detector detected a particle. This 
would absorb a particle from the vacuum, leaving the 
vacuum in a lower energy state, thus extracting energy 
from the vacuum - a form of energy teleportation. In our 
case, the ability to communicate the result to somewhere 
where the partner could be detected would be difficult 
(due to causality) but in some cases [BJ one can actually 
carry out such a procedure and extract energy from the 
vacuum state leaving the system with locally less energy 
than the vacuum (but leaving the system as a whole of 
course with higher energy). 

Let us also look at a more complex example. In this 
case let us assume that the function e(t) has the form of 
a trapezoid - it rises linearly from 0 at time —r to eg 
at time —T, remains constant to time T and then falls 
linearly to zero at time r. The Fourier transform of this 
e(t) is 


state, but it makes, as we saw above, the finding of the 
partner much easier. We will use the zero of {a|/3} (E) 
for E nearest 40, namely E = 38.48966. 



FIG. 1: The evaluation of {a|/3} ( E ) for the trapezoidal cou¬ 
pling as a function of ET with r = 1.2 T. The zeros corre¬ 
spond to orthogonality and pure two mode squeezing. 


For those values of E where {ol\(3} (E) is zero, and 
thus a is orthogonal to /3, the partner mode will be 


«()=/*,[ C, c °s[^ + gFl -cosl^ + Etn 


+c 2 


(w + E ) 2 

cos[(w — E)t\ — cos[(w — E)T] _ it 
(w - E) 2 6 


m 


e(w) 


2eo 


cos(wT) — cos(wt) 
oj 2 {t-T ) 


(36) 


where C\ and C 2 are appropriate normalization factors. 
For example, if we take 


which gives 

f\k\ cos[(|fc| — E)T] — cos[(|fc| — E)t\ 

V 4^ £ ° (|jfe| - E) 2 ’ 

,/W cos[(|fc| + E)T] - cos[(|fc| + E]t) 

V 47t £ ° i\k\ + E) 2 ' (37) 

The requirement that a and (3 be orthogonal can always 
be satisfied for suitable values of E. Their overlap 



{ol\ f3}(E) = J dka* k /3k = 2 J dka* k /3 k , (38) 

o 

is plotted in Figure [l] for r = 1.2 T as a function of ET 
and we see that there are values of ET which make this 
overlap zero. This is certainly not required, as partners 
exist even if one does not have a pure two mode squeezed 


f f cos[(w — E)t\ — cos[(w — E)T]' 

c “ = J (-57V Ef - 1 “ 

= 5.15932, 


Cb = 


/cos[(w + E)t\ — cos[(w + E)T]' 


(w + E) 2 


= 0.0001195, 


(40) 


then, defining tanh 2 r = Cb/c a , we get 


Ci 

c 2 


cosh r 
\[Eb 


sinh r 



= 91.47, 


= 1 . 020 - 10 " 5 . 


(41) 


In this case, the Fourier transform of the partner mode 
does not vanish at w = 0 (because e(E) is not zero), but 
has a step at to = 0. This implies that the partner mode 
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cf>P will have a slow falloff of order l/\t\ for large values 
of t. In Figure [5] we have a plot of the magnitude of the 
partner mode as a function of t for T = 1, r = 1.2 and 
E « 40. In this case, the partner mode is concentrated in 
the same area as is the original detector mode but with 
a far longer tail. However, as we saw above, there is no 
requirement that the partner mode be near the peak in 
the detector mode. 



U 


FIG. 2: The amplitude (f>p as a function of t. 


IV. PARTNERS AND AMPLIFIERS 


An example of a system where the detector or Hawk¬ 
ing mode is completely separate from the partner more 
is the case of a phase insensitive amplifier. Let us take 
the model of such an amplifier as given in Q, in which 
the amplifier is represented as the coupling, by a free sin¬ 
gle degree of freedom q, of two massless one-dimensional 
fields <p and ip, one (ip) having a negative action 


C — Crj) + Cil> + Cq + /lint 

= \ [{M? - (d x cp) 2 ] - i [(d t ip ) 2 - ( d x ip ) 2 ] 


+ ^ [( d t q ) 2 


2q(n<p + uip) S(x). 


(42) 


These fields are supposed to live on the positive x-axis 
with Neumann boundary conditions at the endpoint. In 
order to avoid that the support of the <5(x)-coupling co¬ 
incides with this endpoint, we assume that x £ [—e,oo) 
and consider the limit e 4. 0. The boundary consitions 
then read d x <p(t, —e) = d x ip(t, —e) = 0. 


This model has solutions 


(p(t,x) = cp 0 (t,x) - fiq(t - x), (43) 

ip(t,x) = ipo(t, x) vq(t x) , (44) 

d 2 q + (n 2 - v 2 )d t q = d t [n(po(t,0) + pip o (t,0)\ , (45) 

where (po and ipo are solutions to the free (homogeneous) 
equation (/jl = v = 0) with the above boundary condi¬ 
tions. Note that q is damped as long as v 2 < /r 2 . 

The quantum operators <f>, T, and Q obey the same 
equations. Taking the Fourier transform and expressing 
the solutions in terms of annihilation and creation oper- 


ators, we have 



To(M) = 

oo 

/ -^= (a^e-^ + ate^) 

j v™ 

0 

cos(ocx) ,(46) 

T 0 (t,x) = 

oo 

,/ (‘b-“ + !«") 

0 

cos(wx) ,(47) 

Q(t) = 

que-^ + qte^, 

(48) 

<L = 

2 i /idcj + i 

(49) 

w + i([i 2 - v 2 ) ‘ 


Note that the positions of the creation and annihilation 
operators of the 4' field is reversed since its conjugate 
momentum is = —dtip. As a result, the inner products 
of the two fields [see also Eq. m] 


(4,\i) = 

0 

oo 

(V#) = I Jdx ( ip *n; - 7r^*) , (50) 

o 

are of opposite sign for modes d> u and ip u with the same 
w. The vacuum state for the To field is a maximum of 
the energy, rather than a minimum, and is annihilated 
by the b^ operators. 

The initial (input) fields are those that behave as 
(pin (t + x) or ipin(t + x) while the final (output) fields 
go as (p ou t(t — x) or ip 0 ut(t — x), respectively (remember 
that x > 0). Assuming that the input fields are the free 
fields in Eqs. CTH) and (H71) . 

OO 

T in = T[, n = / -^= fa w e- < " (t+x >+o t e <w(t+ *J') ,(51) 

J V 47 TUJ ' ' 

0 

and similarly for 4h n = 4 /q 11 , the output part is, according 
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to Eqs. (1431) and (l44jl . given by 


^out(t X) 


^out(t x) 


dui 




3 — iuj(t—x) 


7TUJ 


-bt, 


2 i/iv 

uj + — v 2 ) 


duj 


o 

+d u 




e -Mt-x) n 


7TUJ 


UJ — + v 2 ) 

u uj + i{[i 2 — v 2 ) 

+ h.c., (52) 

uj — i(^ 2 + v 2 ) 


2i{iv 


uj — i(/j, 2 — v 2 ) 


uj — i(/j, 2 — v 2 ) 

+ h.c., (53) 


where h.c. is the Hermitian conjugate. Thus, one can 
write the output annihilation and creation operators in 
terms of the input by 


A-u) 

Bp 


. oj - i(n 2 + v 2 ) 
“ uj + i(n 2 — v 2 ) 
f t u-i{n 2 + v 2 ) 
w u - i(^ 2 - v 2 ) 



+ du, 


2 iqLV 

uj + i(fj, 2 — v 2 ) 
2 i[iu 

uj - i(n 2 - v 2 ) 


(54) 

(55) 


Writing the first equation as A u = + /3 w 6j, we see 

that the factor |a w | is larger than unity (unless v = 0) 
which means that signals in the a u channel are amplified. 
However, as required by unitarity, this goes along with 
additional noise stemming from the b J, term. 

Now let us consider a mode in the <f> output channel, 
say, defined by at late times 


f H (t-x)= [ ~^= f H (uj) e _iw ( t_x ), (56) 

J y/Anuj 


where we will assume that fg is non-zero only for uj > 0 
and is normalized so that f duo\fH{uj)\ 2 = 1. The opera¬ 
tor associated with this mode at late times is 


If we define the following frequency dependent mixing 
angles and phases 


cosh du; 
sinh 9 U 


e 




—iuj — (/i 2 + v 2 ) 

—iuj + {n 2 — v 2 ) 

2 fiu 

iuj + (fj , 2 — v 2 ) 

-iuj + (m 2 - v 2 ) 

\ — iu) + (/j , 2 — v 2 )\ ’ 
(. n 2 + v 2 ) + iuJ 
\iuj+(^ + u 2 )\ ’ 


(60) 

(61) 

(62) 

(63) 


then 


sinh d = duj 


f*H («) 


sinh 2 th 


(64) 


Then we find that we can express the annihilation op¬ 
erator of the partner mode in terms of the outgoing cre¬ 
ation and annihilation operators in terms of either the 
input annihilation operators, or of the output 


a p = 


OO 

J duj fA{uj)e~ la ^ (— tanh(i?) cosh(0 w ) 


e lX “ a. 


— coth(i?) sinh(# a 


= J duj [ 2 

o 

sinh(26><J ^ 


(65) 


sinh 2 d — sinh 2 9 u 
sinh(2$) 


sinh(2^) 


( 66 ) 


a-H = = J dufA(u)A u = J dujfp{uj) x 

0 

„ uj — i(fi 2 + v 2 ) - 4 . 2ifu / 

UJ + I^-V 2 ) - “UJ + i(d 2 - V 2 ) 

(We use the' to remind ourselves that this operator may 
be in either of the two channels, before or after the the 
interaction with q) As a result, the partner mode will be 



dp= duj /|f (w) ( du, 


o 


uj — i(fj 2 + v 2 ) 
uj + i{n 2 — v 2 ) 


tanht? 


-bl 


2 ifju 


" uj + i{fj? — v 2 ) 
where the mixing angle is given by 


coth$ 


). 




(58) 


sinh 2 i? = duj 


/■ 




Afj, 2 v 2 


j 2 + {is 2 — v 2 ) 2 


(59) 


If \fjj{u>)\ is a highly peaked function about the fre¬ 
quency uj such that 9 U = d, then the first term will be 
zero and the partner mode, made up entirely of will 
be confined completely to the second output channel the 
i/> channel. However, if k is a broad function (non zero 
over a range of order or larger than fi 2 — u 2 , then the 
partner mode will have support in both the ip and the <f> 
output channels - mostly in the former, but partially in 
the latter as well. 

This will also be true in the Black hole case as well, 
which behaves exactly like this amplifier, with the output 
(j) and ip channels being the modes travelling to infinity 
and those falling into the singularity respectively. For 
highly peaked functions of frequency, the partner is be¬ 
hind the horizon, while for broadly peaked functions of 
frequency, the partner has components both inside and 
outside the horizon. This is another indication of the 
non-linear nature of the partner mode. Since any mode 
is the sum of highly peaked functions, one might expect 
that the a broadly peaked Hawking mode might still have 
a partner entirely behind the horizon, but it does not. 



































V. MOVING MIRROR RADIATION 

Before applying the concept of partner particles to 
what has been taken to be a simple toy model for black 
hole evaporation - the radiation given off by an expo¬ 
nentially accelerated mirror - let us briefly review 

the basic concepts. We consider a massless scalar field in 
1+1 dimensional flat space-time 

a<j) = 0. (67) 

At a point-like mirror with the trajectory x m (f), we im¬ 
pose Dirichlet boundary condition 

<l>(t,x m [t]) = 0. (68) 

In terms of the light-cone coordinates 

u = t — x, v = t + x , (69) 

the general solution of □</> = 0 without the boundary 
condition (1681) can be written as a sum of independent 
left-moving <j)\eh(v) and right-moving 0 r i g ht(n) contribu¬ 
tions <p(u,v) = 4>\ e ft(v) + bright(u). The boundary con¬ 
dition (1681) imposes constraints on these two parts and 
thus the quantum field can be decomposed as 

OO 

/ p —icjv _ p —iu){2r[u\—u) 

dco - ' -a™ + h.c. (70) 

V47TU; 

o 

Here (a™)'*' and denote the initial creation and an¬ 
nihilation operators and the function t[u\ is implicitly 
determined by the mirror trajectory 

t[u] = u + x m (T[u]). (71) 

Hence the mode functions in Eq. G0 automatically sat¬ 
isfy the boundary condition (1551) . For a mirror at rest 
x m = const, we find t[u\ = u + x m and thus these mode 
functions simplify to e~ lLJV — e -M“+2x m ) which just gives 
2 ie~ lut sin(u;[x m — x])e~ tuXni as one would expect. Thus, 
assuming that the mirror is a rest initially, the initial 
vacuum state is determined by 

V,>o CIO) =0. (72) 

Similarly, for a mirror x m = Vt moving with a constant 
velocity V, we get t[u] = it/(l — V). In these cases, no 
particles are created - but with an accelerated motion of 
the mirror (resulting in a non-trivial form of r[u]) one 
can create particles out of the initial vacuum. 

In terms of the light-cone coordinates, the proper ac¬ 
celeration of the mirror x m /(l — x^) 3 ^ 2 can be writ¬ 
ten as ii m /('u m ?) m ) 3 / 2 . Similarly, the red-shift factor 
\/(l + x m )/(l - x m ) simply reads y/v m /u m . Now, if 
we choose the mirror trajectory in such a way that the 
proper acceleration of the mirror is proportional to the 


red-shift factor, an observer at rest sees a stationary ther¬ 
mal spectrum given off by the moving mirror. This situ¬ 
ation corresponds to the mirror trajectory 

g-KUm e -K,(t-X m ) 

t + x m — Vm — — , (73) 

K K 

where k is a proportionality constant which sets the tem¬ 
perature. As a result, the mode functions satisfying the 
boundary condition (1681) are given by 

<j>w(u, v ) = e~ lu>v - exp ji ^ e _K “| . (74) 

In principle, since the proper acceleration of the trajec¬ 
tory (1751) vanishes for very early times t j. — oo, we could 
consider a mirror moving along the world-line (1731) for 
all times. However, to make the initial behavior as sim¬ 
ple as possible, we assume that the mirror is initially at 
rest and starts accelerating along the trajectory (1751) at 
u° = 0 which means = — 1/re, i.e., 


v m (u) 


1 J 1 — nu for u < 0 
K \ e~ KU for u > 0 


(75) 


Consequently, incoming light rays with v < — 1/k are 
reflected by the mirror at rest, i.e., initial waves of 
the form e~ lulv are simply transformed to final waves 
of the form e~ lulu . Incoming light rays in the window 
— 1/re < v < 0 are reflected by the accelerating mir¬ 
ror. In this region, initial waves of the form e~ lulv are 
stretched by the increasing red-shift factor and finally be¬ 
have as exp {iuje~ KU / k}. More generally, an initial wave- 
packet of the form 4>i e h(v) in the region —1/K<w<0is 
transformed to fright ( — e _K “/«), i.e., a final right-moving 
wave-packet in the region u > 0. The remaining light 
rays with v > 0 do not see the mirror at all and thus 
their functional form is unchanged. Hence the null line 
v = 0 is analogous to the black hole horizon. 

Starting in the initial vacuum state (1721) . we can de¬ 
rive the two-point functions in the final state. In order to 
avoid artifacts stemming from the infra-red divergence of 
the massless scalar field in two dimensions, we consider 
the first derivatives of the fields. (This is somewhat simi¬ 
lar to considering the electric and magnetic fields instead 
of the scalar and vector potentials.) As mentioned above, 
the field can be split up into a left-moving ^i e ft(u) and a 
right-moving part fright (u). The correlation between two 
final the left-moving contributions (with vi ,2 > 0) gives 


(0| d v 0ieft(vi)d v 0i e tt(v 2 ) |0> = ^ ^ , (76) 

which just reflects the fact the associated quantum state 
is locally indistinguishable from vacuum (since it has not 
“seen” the mirror at all). 

Considering the correlation between two right-moving 
contributions which have been reflected by the acceler¬ 
ated mirror (with u 1,2 > 0), however, gives 


< 0 | du4 > right('U*i')d u (j) T [gl l t(u2') | 0 ) — 


16-7T sinh 2 («;[ui — u^/2) 


(77) 
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As already suggested by the periodicity in imaginary time 
(KMS condition), this is locally indistinguishable from a 
thermal state with the temperature 

t -~t, (78 > 

The fact that this thermal contribution and the above 
vacuum part are actually just two regions of the same 
pure state results in non-trivial cross-correlations be¬ 
tween these two modes 

(0| ^01eft(ui)5 u ^ r ight(M2) |0) = ^_ KU2 + ( 79 ) 

These results can be generalized to different mirror tra¬ 
jectories v m (u) in a straightforward manner. In this case, 
the mode functions read e~ lLJV — e ~ lujv m {u) an( j 
two-point function is given by 

(0| 4>(ui,V 1 )$(U2,V 2 ) |0) = 

-T Inf, (80) 

47T \[Ul - U m ('U 2 )J[u m ('Ui) - V 2 \) 


VI. PARTNER PARTICLES FOR MIRROR 
THERMAL RADIATION 


Now let us try to determine the partner particles for 
the thermal radiation created by the mirror as an ana¬ 
logue for Hawking radiation. Thus, we define the outgo¬ 
ing Hawking wave function as a linear combination 

of final positive-frequency right-moving plane waves 

OO 

f H (u) = J dn.f H (n)e- mu , (si) 

o 


where //r(fl) is then the Fourier transform of fn(u). Due 
to the restriction to positive final frequencies, the sup¬ 
port of / h(u ) is unbounded, i.e., extends to negative u 
as well. However, for simplicity, we assume that 
lies mostly in the thermal region u > 0, i.e., that 
is exponentially suppressed for u < 0. Alternatively, we 
could consider the case of eternal acceleration of the mir¬ 
ror, where the thermal region extends to negative u. 

The Bogoliubov coefficients can then defined by the 
overlap between these modes m and the initial pos¬ 
itive/negative frequency modes e ±lulv . For the trajec¬ 
tory Cli, these overlap integral can be calculated ana¬ 
lytically in terms of T-functions etc. However, instead of 
using these F-functions, we do the following trick: Ini¬ 
tially, the mode (ISTT) behaved as 

OO 

/&(« < 0 ) = j dtn f H (Q) (-^y n/K , ( 82 ) 

0 

and > 0) = 0. Now, let us consider the following 
linear combinations 


/nW = 


e ±nfl/(2K,) f or v <0 

gT’ t^/(2k) \ KjV \ iQ ‘/ k f or v > o 


(83) 


These linear combinations are chosen such that the func¬ 
tion /q (v) is holomorphic in the entire lower half of the 
complex v plane, i.e., for 3(v) < 0, and has a singularity 
at v = 0 as well as a branch cut from v = 0 to v = oo in 
the upper half. On the other hand, recalling the structure 
of the initial mode functions e ~ luJV , we find that any solu¬ 
tion is exactly decomposed of positive (initial) frequency 
modes if and only if it is holomorphic in entire lower half 
of complex v plane. Thus the linear combination ( v ) in 
Eq. (1831) contains only positive (initial) frequency modes 
for all H, i.e., it corresponds to an initial annihilation 
operator <Zq with a 1 /) |0) =0. Conversely, the other com¬ 
bination J'q (v) is holomorphic in the upper half of the 
complex v plane and thus contains negative (initial) fre¬ 
quencies only, i.e., it corresponds to an initial creation 
operator. Using the symmetry [/^ (v)]* = f^ n (v), we 
find that corresponds to (d™ n )' i \ 

Now, we can decompose the function as a lin¬ 

ear combination of /q(u) which gives 


g+^n/( 2 «)/+(<;) _ g— ttQ/( 2 k) j-( v ) 
-\/2sinh(f2/«:) 


|ku| 


iCl/ k, 


(84) 


for v < 0 and vanishes for v > 0. This enables us to 
directly read off the decomposition of the final Hawking 
mode (|81D into initial creation and annihilation operators 


an = 


OO 

JdClf H (n) an +/3 n (a™ n y , (85) 


with the Bogoliubov coefficients 


an = 


g+7rf2/(2re) 

sinh(fi/K) 


Pn = 


g-7rf2/(2K) 

sinh(f2/«:) 


( 86 ) 


where the denominator \J 2 sinh(f2/re) ensures the correct 
normalization |an| 2 — |/?n| 2 = 1. 

This is the starting point for the derivation of the part¬ 
ner mode dp. Note that the modes fni v ) are orthogonal 
with respect to the usual inner product for the scalar field 


(0i 102) = i 


d.^ PI dp 02 , 


(87) 


where <j>\ dp 02 = PI dp p 2 - p 2 dp p\. For purely right- 
moving modes, we may align the hyper-surface E with a 
null-line of constant u (or even J ~) such that the dE^- 
integral becomes an integration over v and the derivative 

dp simplifies to dv As explained above, the mode func¬ 
tions /q (v) are decomposed of purely positive frequency 
initial waves e~ luv with uj > 0. Conversely, the mode 
functions fni v ) are decomposed of purely negative fre¬ 
quency initial waves e +lwv with uj > 0. As a result, the 
contributions ( v ) and (v) are orthogonal and thus 
the vectors a and (3 are orthogonal, i.e., we have pure 
two-mode squeezing for all ///(H). 
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Using the arguments presented in Sec. IE we find that 
the partner mode reads 


ap= dH/^(fi) \xP fid“ n + x l a u (a(")^ 


( 88 ) 


with the factor \ = aj (3 


X 


2 


Jdn\f H (n)\ 


0 


2 C& 


f cm \f H (H)| 2 /? 


0 


2 

n 


(89) 


If the Hawking mode ///(H) is well localized and peaked 
at a given frequency fi = Ho, then we may approximate 
this factor by \ ~ As a result, the wave function 

of the partner particle fp(v ) contains the following linear 
combination of the modes ( v ) which yields 



e +7rn/(2«)^+^ _ e -nQ/(2 K ) 

-\/2 sinh(f 1/k) 


\nv\~ in,K ,(90) 


FIG. 3: The outgoing detected Hawking mode, its shape when 
traced back to the initial state, and its partner mode. 


for v > 0 and vanishes for v < 0. Thus the wave-function 
of the partner particle is approximately the mirror image 
of the initial form of the Hawking mode (1521) on the other 
side of the horizon at v = 0, i.e., 

OO 

f P (v )« J dnf* H (n) = Fh 

0 

for v > 0 and zero for v < 0. As a result, detecting 
a Hawking particle with, say, fi = O(n) in the thermal 
region at late times ku > 1 yields (up to normaliza¬ 
tion) approximately the same state as creating a partner 
particle with exponentially short wavelengths in the left- 
moving vacuum region at very small but positive values 
of v. Figure [5] is a plot of of a specific outgoing Hawking 
mode (say the mode detected by some detector) and that 
mode in the input state, and its partner mode. 

If we could signal the measurement result of the Hawk¬ 
ing detector at large u > 0 on the right-hand side to this 
vacuum region on the left-hand side, we would (at least in 
principle) be able to extract energy out of this quantum 
state, which is locally indistinguishable from vacuum - 
this is directly related to the concept of “energy telepor¬ 
tation” . However, causality prevents us from signaling 
since these two events are space-like separated. 

Another point is that the cancellation of the contribu¬ 
tions in Eq. (1501) occurs at one frequency fi = Ho only. 
If we consider a small but finite width AH -C Ho, the 
partner mode will also have support in the thermal re¬ 
gion u > 0, i.e., v < 0. For simplicity, let us consider a 
Gaussian wave-packet of the form 


ln(«w) 


(91) 


/ ff (H) = A/ - exp 


(H-H 0 ) 2 


+ jfiw 



(92) 


which is centered around Ho in frequency space and 
around uo in position space. Even though this ///(H) is 
not exactly zero for H < 0, this contribution is exponen¬ 
tially small for AH -C Ho and thus negligible. Similarly, 
the tail of this wave-packet for u < 0 can be made very 
small by assuming itoAfi 1. 

After inserting this form (1921) , the formula for the part¬ 
ner mode fp(v < 0) contains an H-integral whose in¬ 
tegrand vanishes (to lowest order) at H = Ho- Taylor 
expanding this integrand around this zero then yields a 
first-order contribution of the form (H — Ho)///(H) which 
can also be represented by (Afi) 2 [m 0 — dn]///(H). After 
the Fourier transformation (ED, the 8q translates into 
iu and thus the partner wave function acquires a small 
contribution with the same support as the Hawking mode 

f P {u) ~ [ M - u 0 ]f H (u ), (93) 

in addition to the dominant contribution m- Due to 
the term [it — uq] , the two modes are orthogonal as they 
should be. 


VII. CONCLUSIONS 

Perhaps the most surprising conclusion of this paper is 
that the partner particles of the thermal radiation (emit¬ 
ted by a mirror or a black hole) are concentrated in a 
region which is locally indistinguishable from vacuum. 
In the black hole evaporation process, the Hawking par¬ 
ticles emitted at early or intermediate times can be en¬ 
tangled not with some other energetic emission at late 
times, but with final vacuum fluctuations. This weak¬ 
ens the usual argument in black hole evaporation studies 


2(AH) 2 
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which assume unitarity in the above sense, which states 
that the large amount of information left in the black 
hole (entanglement with the emitted thermal emission) 
must be accompanied by the eventual emission of large 
amounts of energy. 

Thus, if one were to imagine a black hole constantly 
fed by a pure state designed to just compensate for the 
energy emitted by the black hole in Hawking thermal 
emission, for 10" times the natural decay lifetime of the 
black hole, there must be a huge amount of information 
inside the black hole, encoded in the entanglement with 
the outgoing Hawking radiation. When the black hole 
eventually evaporates that information, which must be 
emitted at late times, one could expect that it must be 
accompanied by a large amount of energy as well. How¬ 
ever this paper offers the possibility that that eventual 
emission of information could be in the form of the vac¬ 
uum, and carrying no energy. 

For example, the Bardeen model [Ij|, a response to 
the AMPS [l4| argument that either unitarity (as de¬ 
fined above) or the regularity of the space-time at horizon 
must be wrong, has the partner radiation to the Hawk¬ 
ing emission trapped within the apparent horizon of the 
black hole, until eventually that horizon disappears. This 


would seem to require a massive emission of energy just 
at the time when that apparent horizon disappears to ac¬ 
company that massive emission of information (i.e., the 
entangled partner radiation to the earlier Hawking emis¬ 
sion). Our results offer the possibility that those partner 
modes are, in that final stage, simply a part of the vac¬ 
uum state with no energy accompanying them. 
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Appendix A: Uniqueness Proof 

In the following, we show that the ansatz (fT21) for 
the partner particle is the most general ansatz one can 
make, i.e., that the partner particle is uniquely deter¬ 
mined by our two conditions (unless we have pure single¬ 
mode squeezing a||/3). For free fields, the initial vacuum 
state is a Gaussian state and thus the state restricted to 
the two modes dp and dp must also be a Gaussian state. 
In the position representation, i.e., as a function of the 
position vector x = (xp,xp ) where xp = (dp + 4)/v 2 
and Xp = (dp + dp)/\/2, the most general wave function 
of a pure Gaussian state reads 

ip(x) = Af exp | — i x ■ M ■ ;c| , (Al) 

where M is a symmetric but possibly complex matrix 
and Al the corresponding normalization factor. In order 
to have a normalizable state, the real part of M must 
have two positive eigenvalues Ai^- 
Now, the derivative of tp yields 

ip(x) = —M ■ xip(x) . (A2) 

ox 

In terms of the momentum operator p, we get 

(ip + M ■ x) \ip) = 0 . (A3) 

This motivates the introduction of pre-annihilation oper¬ 
ators A = ip + M ■ x which obey the following commu- 



12 


tation relations 


Ai,Aj 

Ai,A\ 


a\,a\\ = 0 , 

Mu + M*j = 25i(Mjj). 


(A4) 


Since 1ft {M) is a real symmetric and positive matrix, we 
may diagonalize it with an orthogonal (rotation) matrix 
D such that D ■ Sft(M) • D * = diag{A/}. As a result, the 
operators a = D [2 Sft ( M )] -1 / 2 • A , i.e., 


ai 


DijAj 

V2A i 


d/ IV 1 ) = 0 


(A5) 


construction above, we see that the Hawking and partner 
mode operators dp and dp must be linear combinations 
of these operators di and <22 as well as their adjoints a\ 
and dj. Thus, the linear sub-space spanned by ai and 02 
can be identified with that of d|| and o_l and we arrive at 
the ansatz m- 

Alternatively, one could insert the general ansatz for 
the modes dp = {a|d} + {a\(3} and dp = { 7 |a} + {a|d} 
into Eq. (I A. 'll) which gives the two linear equations 

(3 - a + Mn(/3 + a.) + M i2 (d + 7 ) = 0, 

5 — 7 + M 2 i(<5 + 7 ) + M 2 2(/3 + a) = 0, (A 6 ) 


satisfy the standard commutation relations and do also 
annihilate the state |?/>). Since this state |?/>) is just the 
initial vacuum state reduced to the two modes dp and 
dp, the two operators di ,2 above must be a linear com¬ 
bination of the initial annihilation operators. ^From the 


where Mu are the components of the symmetric matrix 
M (which also depend on a, /3, 7 , and S). Since these 
two equations are linearly independent for all M/j, we 
find that 7 and 6 must lie in the same sub-space as a 
and (3 (which are assumed to be linearly independent). 









